Can massless QCD dynamically generate heavy quarks? 
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As it was suggested by previous works on a modified perturbation expansion for QCD, the possibility 
for the generation of large quark condensates in the massless version of the theory is explored. For 
this purpose, it is firstly presented a way of well define the Feynman diagrams at any number of 
loops by just employing dimensional regularization. After that, the calculated zero and one loop 
corrections to the effective potential indicate a strong instability of the system under the generation 
of quark condensates. The also evaluated quark condensate dependence of particular two loop terms 
does not modify the instability picture arising at one loop. The results suggest a possible mechanism 
for a sort of Top Condensate Model to be a dynamically fixed effective action for massless QCD. 
The inability of lattice calculations in detecting this possibility could be related with the limitations 
in treating the fermion determinants. 



I. INTRODUCTION 

The explanation of chiral and flavor properties of QCD 
is one of the fundamental research issues in high energy 
physics P, 0, 0, 0, IS El • A central problem in this field is 
the question about the origin of the quark masses, which 
indirectly determine the nature of most of the observable 
matter. More generally, the explanation of the elemen- 
tary particle mass spectrum is considered as a fundamen- 
tal question posed on the research in Physics. This situ- 
ation gives relevance to the close examination of possible 
mechanisms that could be playing a role in this problem. 

Specifically, a modified perturbative QCD, altered in 
a way that incorporates the presence of a condensate of 
zero-momenta gluons and quarks in the initial state used 
for constructing the Wick ex pan sion, have been consid- 
ered in Refs. S H ES El lH E An interesting as- 
pect of these works was the BCS-likc modification of the 
gluon state studied in Refs. Q, H, which directly led to 
the prediction of the constituent masses for light quarks, 
after simply fixing to its currently estimated value, a rel- 
atively unrelated quantity, the gluon condensate parame- 
ter. Therefore, it looked reasonable to expect that a sim- 
ilarly constructed state for the quarks, after to be intro- 
duced in the treatment of massless QCD, could have the 
chance of generating the Lagrangian quark masses. The 
result of the preliminary consideration of this idea in Rcf. 
ITS was positive. A modification of the free-quark prop- 
agator by introducing the zero momentum terms repre- 
senting analogs of the Cooper pair condensates within 
this problem was done. After that, the simplest approxi- 
mation for the Dyson equation for quarks produced a di- 
agonal Lagrangian mass matrix by the choice of an also 
diagonal structure, for the quark condensates. Hence- 
forth, the results of the work E, pointed out that the 
presence of colorless quark anti-quark condensates in the 
vacuum of the free theory generating the Wick expan- 
sion, is able to produce the observed quark mass matrix. 
Therefore, the question arose about the possibility for 



the generation of the necessary pattern of quark conden- 
sates by a dynamical breaking of the flavor and chiral 
symmetries. A preliminary step in the study of this issue 
was also done in Ref. |13| . where the one loop contribu- 
tion to the Cornwall- Jackiw-Tomboulis (CJT) effective 
potential for composite operators EH was evaluated as 
a function of the parameters determining the quark and 
gluon condensates. However, it should be underlined that 
as the CJT potential is not dir ectly giving the ground 
state energy at its minimum |l4l Il5| , the implications of 
the results of the evaluation for the dynamical symmetry 
breaking problem under study were not clear. We esti- 
mate that in the case that a similar mechanism to the 
one acting in standard superconductivity, could be play- 
ing a role in the problem, the generation of large quark 
condensates could be expected to be produced by strong 
binding color forces, linked with the interaction vertices. 
On another hand, the inability of the standard perturba- 
tion expansion in evidencing this effect, looks for us to be 
rooted in the explicit disregarding of the inclusion zero 
momentum quarks and gluons in the free vacuum state, 
before the adiabatic connection of the interaction. The 
incorporation of these modified Lorentz invariant vacu- 
ums in a BCS style, allows to produce non-trivial modi- 
fications in the Wick expansion as it was argued in Ref. 



In the former works on the theme |a|alalifllliUJ,llii 

a troublesome aspect was remaining about the appear- 
ance of singularities in the Feynman diagrams due to the 
presence of Dirac Delta terms in the free propagators. 
This issue will be approached here in one of the Appen- 
dices. The basic idea of the adopted procedure is to sim- 
ply be consistent with the dimensional regularization and 
to extend the appearing (5(0) like singularities for continu- 
ous D dimensions in the way early introduced by Capper 
and Liebbrandt [l^j . This procedure lead to the outcome 
that these factors simply vanish in the D — > 4 limit. An 
additional recipe is also taken for the evaluation of the 
l/(p 2 + ie) factors at zero momentum. These terms also 
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appear due to the joining of more than n — 1 different 
condensate lines in a vertex having n legs. The rule cho- 
sen in this case will be to just regularize the scalar field 
propagator at the zero value of the 4-momentum point 
to be equal zero. These points are discussed in Appendix 
A. The selected prescriptions, in addition, allow to iden- 
tify the propagators evaluated in |l3j |. as modified tree 
propagators. They also have zero order in the coupling 
g series expansion after considering as the independent 
parameters of the theory, the proper <?, in addition to the 
gluon and quark condensate parameters multiplied by g 2 . 
This transformation seemingly will allow to rearrange the 
full loop expansion to a form in which the propagators 
derived in Ref. will play the role of new tree propaga- 
tors. However, this more formal aspects will be relegated 
for a further study. Here we only consider the one and 
two loop cases in which it is clear that the summation 
over the zero order in the coupling self-energy insertions 
can be done. The full demonstration of the coincidence of 
the terms evaluated here with the exact loop expansion 
terms of the alternative series needs only for the checking 
about wether the combinatorial factors of all the lower 
loop terms arising from the higher loop ones could ob- 
stacle the proof. The verification of this issue will be 
considered elsewhere. 

After giving the procedure for defining the perturba- 
tive expansion, we continue in this work the study of 
the possibility for the dynamical generation of masses 
in massless QCD. Zero and one loop vacuum contribu- 
tions to the effective potential are evaluated and also the 
full dependence on the quark condensate of a particular 
two loop contribution is presented. As noticed above, 
the usual definition of the effective potential was evalu- 
ated here (|l5j) m order to avoid the non direct bounded 
from below property of the directly evaluated CJT po- 
tential p^ |. The approximation considered consists in 
inserting all the condensate dependent parts of the one 
loop self-energy corrections into the free propagators, in 
the usual zero and one loop correction, as follows from 
the discussion in Appendix A. This procedure results in 
employing in the usual loop diagrams, the propagators 
that produced the constituent masses for light fermions 
in the work Q, by also considering propagators associ- 
ated to the gluon and quark condensates. The quark 
condensate dependence of only one particular two loop 
diagram, was also evaluated to estimate the possible ef- 
fects of the next corrections. The evaluation of the full 
two loop dependence on the condensates, in order to de- 
termine their net effect on the results will be considered 
in further studies. We prefer to postpone the calcula- 
tion of these terms, waiting for a precise definition of the 
renormalization scheme in the modified expansion. 

The results obtained here for the effective potential 
indicate a dynamical generation of quark and gluon 
condensates. The dependence on the potential on the 
quark one, get an unbounded from below behavior, in 
the present approximation. This instability becomes 
stronger by increasing the gluon condensate. This point 



suggests the relevance of the presence of the gluon con- 
densate for the dynamical generation of the quark one. 
The unbounded from below dependence, then indicates 
the need of higher approximations for producing an even- 
tual minimum of the potential. The picture arising pre- 
cisely reproduce the one expected to occur in Ref. [l3| 
as a possible consequence of the underlined analogies be- 
tween the construction of the modified wave functions for 

QCD and the BCS states HEl. The de P endence of the 
potential on the gluon condensate starting at zero value, 
is a decreasing one, which become steeper when the quark 
condensate grows. Therefore, this property generate the 
expectation about that, the yet to be determined stabiliz- 
ing contribution to the potential naturally could show op- 
posite behavior in the quark and gluon dependence, able 
to produce not only the stabilization, but also a value for 
M at the minimum laying near the known constituent 
mass M = 333 MeV. At low values of the quark con- 
densate parameter, the gluon one develops a minimum. 
However, in the present approximation, for fixed values 
of the quark condensate the value of the gluon conden- 
sate at the minimum tend to grow, when the value of 
X increases. This is not unexpected, since we have not 
evaluated the full two loop dependence on the gluon con- 
densate, nor the possible existing stabilizing terms. They 
as mentioned before, could help in maintaining the gluon 
condensate at a minimum being near low values, allowed 
to be fixed to the observable value. This fixation could 
be done, let say by selecting scale parameter ji. It will be 
clearly surprising that the magnitude of the X required 
for the stabilization could be as high as of the order of 
hundreds, as it would be needed for predicting the top 
condensate mass near mtop = X = 175 GeV (as given by 
the pole of the quark propagator for high X values [l3jp. 
However, we could not yet disregard this possibility and 
the search for an estimation of the stabilizing terms will 
be undertaken. 

A complementary dependence of the two loop potential 
as a function of the quark condensate has been evaluated 
in the form of a simple 2D integral depending on this 
condensate and the gluon condensate one. The outcome 
for this quantity turned to be finite after including the 
corresponding quark condensate dependent part of the 
usual quark counterterm. This result gave us confidence 
about that the renormalization procedure can be well 
implemented in the modified theory. However, a careful 
discussion of this question should be considered in detail. 

The introduction of the gauge parameter dependence is 
another problem which need an additional careful consid- 
eration. Partial argues about the gauge invariance of the 
scheme have been done 0, 0] , and also the evaluation 
of the one loop gluon self-energy directly satisfied the 
transversality Ward identity 9J. The discussion in the 
present work and in Ref. [IJ make also clear that the 
modified theory is a multi-parameter one in which the 
implementation of the gauge invariance can show sub- 
tleties needing to be carefully addressed. However in the 
appendices we present some ideas about how consider 
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this questions. We also expect to be able of presenting 
more definitive conclusions in next studies. 

The work is organized as follows: in Section 2, the 
propagators employed in the calculations are presented 
and the effective action vacuum diagrams described. The 
zero and one-loop potential contributions are discussed 
in Section 3. Section 4 consider the quark condensate 
dependent part of the two loop effective action contri- 
bution being considered. The Section 5 is devoted to 
expose and discuss the results of the evaluations done. 
Appendix A, exposes the procedure for eliminating the 
singularities in the diagram expansion through employ- 
ment of dimensional rcgularization and Appendix B dis- 
cuss gauge invariance aspects and the perspectives of its 
implementation in the proposed scheme. In the summary 



the main results of the work are shortly reviewed 



II. PROPAGATORS AND EFFECTIVE ACTION 

In next sections the evaluation of the effective potential 
including zero, one and two loop corrections will be con- 
sidered. The contributions will be calculated by insert- 
ing the infinite ladder of condensate dependent one-loop 
self-energy parts in the original free propagators follow- 
ing the rules defined in Appendix A and Ref. ^3| . These 
propagators for quarks and gluons, as well as for the con- 
densate lines (defined in Appendix A and Ref. 0|) are 
given as 
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where l|ll5|l are the gluon, quark and ghost propagators 
respectively and (0} , (J2J) the gluon and quark condensate 
ones. In this work we will adopt the general conventions 
for the spinor, color and Lorentz groups, the free propa- 
gators and interactions vertices of the reference [18| . 

The parameters to 2 , M, Sf are related with the con- 
stants C and Cf (See Ref. characterizing the gluon 
and quark condensates, as follows 
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In these relations the parameter / = ±(§) 2 will be 
considered for the two values of its sign. However, the 
negative value was the only one implied by the calcu- 



lation done in 13] . In that work, it followed that if 
the parameter C is chosen as positive, then the gluon 
mass in the simplest approximation turns to be tachy- 
onic (to 2 = — ) and on another hand the constituent 
mass for light quarks M becomes real. The situation re- 
verses for negative value of the gluon parameter C. How- 
ever, it should be taken into account that these relations 
appeared due to the special fact that the condensate de- 
pendent part of the one gluon self-energy had no contri- 
bution from the fermions. This situation occurred pre- 
cisely because they were assumed as massless. Therefore, 
as the very same motivation of the present discussion is 
related with the possibility of generating masses for the 
quarks, and also paying attention to the fact that the cor- 
rected propagators could be perhaps also constructed in 
a self-consistent manner (allowing for a self-consistently 
generated mass) , we consider of interest to also examine 
the calculation for the positive choice of / . It should be 
mentioned that for this value of the parameter, both the 
gluon mass and the constituent masses are real. More- 
over, the gluon mass is near the value estimated though 




FIG. 1: The figure shows the seven Feynman diagrams defin- 
ing the zero and one loop contributions evaluated in the work. 
The lines having cuts correspond to the condensate propaga- 
tors. Therefore, although the associated diagrams may look 
as two loop ones the Delta functions associated to them ef- 
fectively cancel one of the loop integrals. 



other studies ^3] m = 0-5 GeV, while the constituent 
mass gets the reasonable value M = 0.33 GeVp|. How- 
ever, a study of the physical justification of the positive 
choice of / needs to be done. This question could be 
considered elsewhere. 

It should be explicitly stated that only one flavor was 
assumed to be condensed en the present discussion. This 
condition was elected because at present level of approx- 
imation , the consideration of various flavors will simply 
lead to an of addition identical fermions contributions 
to the potential. The question about the possible inter- 
ference of various quark condensates, since it needs for 
higher order approximations for its appearance, will be 
relegated to further studies. It is clear that this is a rel- 
evant point, because only in the case that the presence 
of various kinds of such condensates will be rejected by 
the system the dynamical generation of only one (main) 
quark condensate will be preferred . This could occurs by 
example, due to the presence of terms in the effective po- 
tential growing in value when more than one condensate 
are present. This effect seems clearly possible to occur 
but its considerations need for at least three loop correc- 
tions (or their descendants according to the reasoning in 
Appendix A) in which different quark loops can start to 
appear |l3[ . 




FIG. 2: The figure illustrates: a) The two fermion contri- 
bution Fa considered in the work for getting a sense of the 
influence of the higher loops; b) The substracted diagram 
which is the same T% evaluated at S — 0; c) The diagram 
associated to the fermion counterterm Tio and d) This last 
diagram taken at S — indicated by fu. 



The collection of zero, and one loop diagrams which 
were evaluated are illustrated in FigQ The diagram Ti 
is the only non having closed loops, that is a tree cor- 
rection. T2 and T3 show the usual gluon and quark one 
loop corrections associated to the propagators Q and 
© respectively. Further, diagrams T4, r 5 ,r 6 and are 
related with the one loop corrections being "descendant" 
from the two loop ones due to the cancellation of one of 
the two loop integrals by a condensate propagator, and 
the insertion of all coupling g independent self-energy in- 
sertions leading to the propagators ^ ID m the other two 
lines (See Appendix A). 

Finally in Fig. [21 the diagrams Tg and Fg, Tio and 
Tn are defined as follows: Tg is the standard diagram 
for the two loop correction including all the coupling in- 
dependent self-energy insertions in its internal lines, Tg 
is the same contribution as Tg but taken in the limit 
S — > which is substracted in order to consider only the 
quark condensate dependent part of the two loop term, 
the limit S — > is indicated by the rings in the fermion 
lines. Finally Tio is the g 2 contribution associated to 
the fermion counterterm and Tn is the same contribu- 
tion in the limit S — > substracted in order to again only 
consider the quark condensate dependent part of the po- 
tential. The momentum integrations, in writing the di- 
agram expressions, will be taken in Minkowski space for 
after perform the Wick rotation. However, it should be 
made precise that in order to make the rotation with- 
out encountering poles, the sign of m 2 should be positive 
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(/ = (f) 2 )- However, we will perform the Wick rota- 
tion for the two signs of / without including the terms 
that could be incorporated by rounding the poles in the 
p variables when deforming the integration contour if / 
takes its negative value. Therefore the results obtained 
for / negative, should be interpreted as the evaluation of 
the effective potential in Euclidean field theory. That is, 
the evaluated quantity corresponds the thermodynamical 
effective potential in the limit of zero temperature. 

The employed expression for the fermion renormaliza- 
tion constant is given by 
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As it was remarked before, only one quark flavor will 
be considered for the present qualitative discussion, since 
up to this level all the quark flavors will produce a sum 
of contributions all of the same functional dependence on 
their respective condensates. However, the fact that light 
quarks exist furnish a guiding principle in the sense that 
the mass acquired by them, if their quark condensates do 
not develop, should coincide with the parameter M 
at the minimum of the effective potential. 

The value of the renormalized gauge parameter [3 is 
equal to one in accordance with the fact that the modi- 
fication of the Feynman rules induced by the presence of 
the condensates was obtained in [8| within the Feynman 
gauge (3 — 1. As it was mentioned in the Introduction, the 
question of the gauge parameter invariant formulation 
of the modified perturbation theory under study, which 
clearly represent a required step in the formal completion 
of the approach should be further considered. However, 
some general remarks on this question are given in Ap- 



pendix B and a more concrete study is expected to be 
considered elsewhere. 



III. ZERO AND ONE-LOOP TERMS 

In this section the results for the evaluation of contri- 
bution to the one loop effective potential T\ to will 
exposed below in consecutive order. 



A. Zero loop term 

The direct substitution of the gluon condensate prop- 
agator m the analytic expression associated to Ti, 
after evaluating all the Lorentz, spinor and color traces 
leads to 



r(o) = 



2m 4 
1? 



.y(0). 



That is, a positive potential proportional to m 4 . As the 
one loop terms have zero order in the coupling g, in the 
expansion in powers of the parameters defined in Ap- 
pendix A, this term shows a power -2 of g, since the 
original diagram was of order two and there is two con- 
densate lines in the diagram which reduce the power in 
four according to (|A3|) . 



B. Standard one loop terms 

The sum of the one loop terms corresponding to T 2 
and T3 in FigQJhave the form 
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which have been expressed as the sum of a S independent 
term corresponding to the same diagrams evaluated at 
5 = plus a S dependent contribution vanishing in the 
limit S — > 0. After calculating the Lorentz, spinor and 
color traces for the S = gluon and quark loops Fg 
and , and dimensionally regularizing the integral, it 



follows 
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But, in both cases the recourse of taking the derivative 
of the expressions over the parameters in the gluon and 
quark cases leads to simpler expressions. Then, after also 
performing the Wick rotation in the temporal momentum 
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component according to 



Pa 



l p 4 , 



the derivative over the parameters expressions can be in- 
tegrated in momentum space by employing the formula 
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in which L can be selected as m or M as appropriate 
for the gluon or quark terms respectively. The results of 
the integrals after to be integrated over the parameter 
again from their zero values to the original ones, by also 
considering that the extended dimension D is such the 
real part of D — 2 is positive, allows to obtain 
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At this point after considering relations ©, © and 
© defining m and M as functions of the dimension 
D , and substracting the pole part in e of i|10|) and I jllj l. 
the Minimal Substraction result for the one loop effective 
action is 
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The S dependent correction after all the trace evaluations can be written as 
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This integral after the Wick rotation is convergent in the limit D — > 4 and takes the form 
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C. One loop terms descending from the two loop gluon diagrams 

After writing the analytical expressions for the diagram T4 and evaluating the Lorentz, spinor and color traces, the 
expression can be rewrote in the form 
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In the second line of this equation the Wick rotation has been done. The integrals can be explicitly performed to 
give 
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After applying the same procedure for the analytical expressions associated to the diagram the result is 
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It is an interesting outcome that after adding these two 
contributions and removing the dimensional regulariza- 
tion limit, the result remains finite, a fact that also delete 
the logarithmic terms in the outcome. The total contri- 
bution of these terms for the potential at the end takes 



the form 



lim (V, 
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D. One loop terms descending from the two- loop quark diagram 

The last one loop diagrams Tg and Tj correspond to the descendants of the two loop terms having a closed fermion 
line. The integral expression obtained for them after performing the Lorentz, spinor and color traces are not so simple 
and we just numerically evaluate them in this work. The resulting integral expressions are 
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In ending this section it can be remarked the interesting outcome that all the "descendant" diagrams became finite 
ones. I 



IV. TWO-LOOP QUARK TERM 



As we are interested in the studying the generation of 
the quark condensates we will consider here only the full 
dependence on this quantity associated to the diagram T 8 
in Fig. [3 defined by a fermion loop formed with propa- 
gators Pjl. showing two quark-gluon interaction vertices. 
Therefore, the same diagram expression but evaluated at 
5 = will be substracted from this contribution. This 
terms is associated with Tg in Fig. [3 This substrac- 
tion simply corresponds to the same analytic expression 
of the diagram but taken for S — and is represented 
by the same figure but having small rings in the quark 
propagators. As the diagram associated to the fermion 
counterterm of the standard massless QCD Fio(of order 



g 2 and therefore needed for renormalization at one loop 
level) is also depending on the condensate parameter S, 
the same kind of substraction is done of the 5 = coun- 
terterm term associated to Tn, in which again the ring 
in the quark line means the evaluation in 5 = 0. 

The substracted terms, exactly give the full two loop 
term formed by two quark propagators and one gluon line 
of the theory in the absence of the fermion condensate. 
As mentioned before we will relegate the evaluation of 
the full two loop gluon parameter dependence to further 
studies. The main reason for doing so is that for these 
terms, it is more relevant to precisely define the way in 
which the renormalization should be done within the con- 
sidered scheme. As it is also following in Appendix A, 
at the two loop level there will appear additional quark 



8 



condensate dependence coming from two loop diagrams 
being descendant from higher loop terms of the original 
expansion as discussed in Appendix A. From exploring 
evaluations we know however, that it seems possible to 
cancel the two loop infinities by renormalizing the con- 
densate parameters. However, a clearer understanding 
on the structure of the allowed counterterms in the mod- 



ified expansion is yet desirable before evaluating the full 
two loop terms. 

After evaluating the spinor and color traces in the an- 
alytic expressions corresponding to the Feynman graphs 
appearing in Fig^ the considered contributions to the 
effective potential can be written in the form 



J 



rg(M,5) 



rg(M,o) 



rg(M,o) 



(A 2 - l)g 2 



d u q d D q' 



1 



{2n) D i (2n) D i {{q - q>) 2 - m 2 ) 
1 



{q 2 (q 2 - M 2 ) 2 - S 2 ){q' 2 {q' 2 - M 2 ) 2 - S 2 ) 
{ - 4q 2 q' 2 (q 2 - M 2 )(q' 2 - M 2 )[(D - 2)q.q> - (3— 



(?' - q) 2 
(N 2 -l)g 2 f d D q d D q' 



W - qf 



(-4) 



(15) 



{q> - q) 2 



(M , _ w- g)g ^- g) )] + 4{D _ w q2q% 



(2ir) D i (2n) D i {{q - q') 2 - m 2 )(q 2 - M 2 ){q' 2 - M 2 ) 



[(D-2)q.qi ~[3 — 
AN{Z 2 - 1) 



M-q 1 



2q.{q' - q)q' .{q' - q) , 



(q'-q) 2 ^ (q'-q) 2 
d D q (q 2 ) 2 {q 2 -M 2 ) 



= 4N(Z 2 -1) 



(2w) D iq 2 (q 2 - M 2 ) 2 - S* 2 ' 



d u q 



(2Tr) D i {q 2 -M 2 )' 



(16) 

(17) 
(18) 



r 



It can be noticed that mass dimension of the parameter 
S is equal to 3, that is a relatively high value. Therefore, 
the terms of the expansion in powers of S for the de- 

(2) 

nominator of the integrand associated to T J will have 



J 



three powers of the momentum convergence factors for 
each power of S appearing in the expansion. The same 

(2) 

effect is occurring in the fermion counterterm 
Then, it follows that the quantity 



T fg (m,M, S,e)=T$(M, S,e) 



ri$(M,0,e) 



T%(M,S,e) 



rg(M,0,e), 



r 



which contains, by construction, the whole dependence of 
the effective action on the fermion condensate parameter 
S, turns to be finite in the limit D — > 4 (e — ► 0). This 
result is simply expressing the fact that the renormal- 
ization constant Z 2 of the massless QCD (determined in 
the absence of any condensate) is also able to extract the 
infinities from the single fermion condensate dependent 
contribution under study. As noticed before, according 
to the above described substraction procedure, the sub- 



stracted terms in addition with the non considered two 
loop ones, exactly correspond to the two loop plus coun- 
terterm contributions in the absence of the fermion con- 
densate. This terms, including the ones descending form 
the higher loops (according to the reasons given in Ap- 
pendix A) will not be considered here. 

The finite contribution T f g before passing to Euclidean 
variables can be written as the sum of the following three 
terms 
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r (2) _. 

L fg ~ 

-(2,1) . 

fg ~ 



r (2,2) 

1 fg 



-(2,3) 



fg 



r (2,l) , r (2,2) r (2,3) 



/9 ^ x /ff 

(TV 2 - l).g 2 



(2tt) d i 



d D q' 



4 (£> - 2) g.g' 



■)x 



(2tt) I5 z ((g - g') 2 - ™ 2 )(<?' 2 - M 2 ) 
S 2 



(4tt) 2 e ' {q 2 (q 2 - M 2 ) 2 - S 2 )(q 2 - M 2 ) 



+ 



(N 2 l).g 5 



d°q d D q' 



2 (D — 2)q.q'S 4 



{2n) D i {2ir) D i (q 2 (q 2 - M 2 ) 2 - S 2 ) 
1 



(g' 2 (g' 2 - M 2 ) 2 - S 2 ){q 2 - M 2 )(q' 2 - M 2 )((q - q') 2 - m 2 ) ' 
P m 2 {N 2 - l)g 2 



d u q d u q' 4 q 2 q' 2 (q 2 -M 2 ) 

X ~ — ^r~, — ~ ~l — rrrm ~t~" X 



4 J (2ir) D i (2ir) D i (q 2 {q 2 - M 2 ) 2 - S 2 ) 

{q' 2 -M 2 ){2 q 2 q' 2 -q.q'{q 2 +q 12 )) 



(g' 2 (g' 2 - M 2 ) 2 - S 2 )((q - q') 2 - m 2 )((q' - q) 2 ) 2 



pm 2 {N 2 -l)g 2 f d D q d D q 



f d D q d D q' 
J {2n) D i (2ir) D -, 



4 g 2 g' 2 S 2 



) D i (g 2 (g 2 - M 2 ) 2 - S 2 ) 



S 2 



(g' 2 (g' 2 - M 2 ) 2 - S 2 ){{q - q') 2 - m 2 ){q' - q) 2 ' 

I 



(19) 
(20) 



(21) 



(22) 



where the term showing the | factor is associated to the 
fermion counterterm. It is responsible for the substrac- 
tion of the divergent part of the remaining expressions. 

After performing the Wick rotation, it is possible to 
eliminate the pole term in e by using the identity 



1 



( 47 r)fr(f)H 2 ] : 



d D q' 



1 



e e£(f,2-f)B(§ - 1, § - 1) J (2ir) D i (q - q') 2 q' 2 ' 

Then, the finite fermion condensate dependent contribu- 
tion to the particular two loop term evaluated here, in 



the limit e — > 0, can be expressed as follows 



Vfg 

vo 



ffg 



-v 



vf+vf+vf+vf+vf 



■(N 2 - l)g 2 M 4 , 



The quantities v 



1, 2, 3,4,5 appearing above were 



reduced to simple 2D integrals after performing the an- 
gular integrations in the 4-dimensional Euclidean space. 
They take the explicit forms 



v? = 2 X' 



POO POO 

• dq dq' 
Jo Jo 



<zV 4 



(g 2 (g 2 + l) 2 + X^)(q' 2 (q' 2 + l) 2 + X*) 
q 2 + q' 2 + 2qq' + f - ie' 



In 



q 2 + q' 2 - 2qq' + f - ie 



,( 2 ) 



+ X' 



POO pOO 

/ dq dq' 

JO Jo 



q 3 q ,3 



1 



(g 2 (g 2 + l) 2 +X B )(g 2 + l) l q' 2 (q' 2 + l) 
q 2 + q' 2 + f-ie , / q 2 + q' 2 + 2qq' + f - ie 



+ 



4 qq'(q' 2 + l) 

q 2 + q' 2 - ie 

a n ln 

4 qq' A 



X ln 



q 2 + q' 2 - 2qq' + f - ie 



q 2 + q' 2 + 2qq' - ie 
q 2 + q' 2 — 2qq' — ie 



}■ 



(3) 



= -X" 



{-1+- 



dq I dq 



g 3 g' 3 



(q 2 {q 2 + l) 2 + X e ){q' 2 (q' 2 + l) 2 + X e )(q 2 + l)(q' 2 + 1) 



q' 2 — ie 



4 qq' 



ln 



q 2 + q' 2 + 2qq' - f 2 - ie 
q 2 + q' 2 - 2qq' - f 2 - ie 
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,( 4 ) - 



dq / dq' q 3 q' 3 (q 2 + l)(q' 2 + 1) x 
Jo 

q 2 g ,2 



{ 



{q 2 {q 2 + l) 2 + X e )(q' 2 {q' 2 + l) 2 + X 6 ) (q 2 + l) 2 (g' 2 + l) 2 



{1 



{q 2 + q 12 



(g 2 -g' 2 ) 2 
f 



4 ?g' 



■In 



(g - g') 2 + / - (g + g') 2 - ^ 

(g + g') 2 + / - (g - g') 2 - «e 



.,(5) 



(3X e 



dq / e?(/ 



4 /4 
g g 



(<7 2 (<Z 2 + l) 2 + XS)(q' 2 (q> 2 + l) 2 + A 6 ) 
(g-g') 2 + /-^ (g + g') 2 -*^ 



In 



(g + g') 2 + / - ie (<? - g') 2 - * e 



r 



in which as before the dimensionless quantities X is given 
as follows 

5 = M 3 X 3 . 

The e parameter is retained here since it helps to regu- 
larize the integrals even in the Euclidean case when / is 
negative. 



V. DISCUSSION 

In this section we will present the results for the eval- 
uation of the effective potential as a function of the con- 
densate parameters M, S and the couplings constant 
g. The calculations were done for the two signs of the 
parameter / that defines the relation between the con- 
stituent quark and gluon mass parameters m and M 
through 

m 2 = / M 2 . 

As it was remarked before, only the negative sign was 
arising in the work because the constituent mass value 
evaluated in that work was satisfying the above relation 
with the negative sign. This fact was a direct conse- 
quence of the free standard quark propagator being mass- 
less. However, as it was noticed before here, we suspect 
that a sort of self-consistent treatment could lead to an 
unrestricted sign of / . Thus the evaluation for positive 
/ values was also considered. An interesting point in 
this sense, is that for positive / the results for the poten- 
tial are completely real, a fact that is not occurring for 
the more relevant case under study that is / = — (|) 2 . 
Moreover, in this situation, it turned out that the value 
of m following, once the quark condensate < g 2 G 2 > is 
fixed, is m = 0.5 GeV , which coincides with a result ob- 
tained in Ref. 17]. Unfortunately, the sign of the light 
quark masses which follows from the Dyson equation 
is opposite to the sign of the gluon condensate parame- 
ter, that is basically the sign of m 2 . Therefore, in case 
that we select the negative sign of / , then the absolute 
value of the constituent mass for light quarks will be also 




FIG. 3: The landscape picturing the dependence of the func- 
tion V (the effective potential divided by l/(87r 4 )) on the 
quark condensate (measured by X) and the gluon one (mea- 
sured by the constituent mass M). The values of the other 
parameters are g — 2.74 and /i = 6.8 GeV. At X = the po- 
tential develops a minimum at certain value of M which can 
be varied by changing the scale parameter fi. The dependence 
on X indicates an instability upon the generation of values of 
X which is not controlled at large X values. 



333 MeV but the mass will be tachyonic. This perhaps 
is another possibility which could be needed to be also 
examined. In any case, neither gluons or quarks appear 
in Nature and perhaps both will be absent as real excita- 
tions in both descriptions in which none of them will be 
asymptotic states after including more corrections [l3|. 

Let us define for the graphical illustrations the 
quantities V(X, M, g, ji) and its imaginary part 
Vi m (X, M,g, n), as the sum of all the contributions to 
the effective potential (the negatives of the effective ac- 
tion terms) evaluated in previous sections divided by the 
constant factor l/(87r 4 ). In the various figures below the 
dependence of the quantity V or its imaginary part Vi m 
are plotted as functions of two of their arguments select- 
ing the others as given by characteristic values of interest 
in the present state of the discussion. The plots are as- 
sociated to the relevant case / = — (3/2) 2 and sometimes 
comments about the effect of the graphs of changing the 
sign of / will be done. The Fig. [3] illustrates the be- 
havior of the effective potential as a function X and the 
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FIG. 4: The dependence of the function V on X and M for 
a wider range of values of X and the same fixed parameters 
g = 2.74 and /i = 6.8 GeV. It shows more explicitly that the 
system in the taken approximation does not show a minimum 
of the potential when X grows. It can be also seen that even 
at zero value of the gluon condensate (M = 0) the potential 
remains unbounded from below. That is, the sole presence of 
the quark condensate also makes the system unstable under 
the generation of X from the state at X — 0. This property 
suggests that the instability effect could not be destroyed by 
the deconfinement transition at high temperaures. 



constituent mass M. both quantities are defined in the 
text in terms of the quark condensate and the gluon one 
through 



St 



_ = 1 ( 9Cf c,) 



X 

M = , for M real. 



/ 

The value of the coupling g selected for the plot was 
g = 2.74 which corresponds to a strong coupling value 

2 

a = j- being near 0.6. In addition the mass scale pa- 
rameter value fj, = 6.8 GeV was fixed. Note that the 
minimum at zero quark condensate X = is laying near 
200 MeV, which is lower but near the constituent mass 
value M = 333 MeV. It is interesting that to fix min- 
imum of the potential for X = 0, at this value of M 
requires a relatively large value of /i. 

As it can be observed, the landscape of the potential 
makes clear that the system at X — dynamically de- 
velops a gluon condensate parameter with a potential 
similar in form to the Savvidi one in the early Chromo- 
magnetic field models. (113, S3)- ^ also can be seen that 
the system at zero values of both parameters shows an 
instability upon the generation of both gluon and quark 
condensate. The instability is stronger for the dynami- 
cal generation of the quark condensate. It can be also 
observed that the increasing of the gluon condensate pa- 



FIG. 5: The plot of the function Vim for the same range 
of X and M used in Fig. 2] and the same fixed parameters 
g = 2.74 and ju = 6.8 GeV and / = -(3/2) 2 . It can be seen 
that the ratio between the imaginary and the real part of the 
potential decreases for the higher values of X. For / = (3/2) 2 
the potential is real. 



rameter makes stronger the instability to the generation 
of the parameter X. This property is supporting the ex- 
pectation expressed in [3, [l3 about that the color cou- 
pling could produce a sort superconductivity effect being 
able to generate intensive quark condensate values alike 
to the Ginzburg-Landau fields. If such effect is really 
occurring in Nature, the Top Condensate model could 
emerge as a possible effective field theory determined by 
the strong forces and upon this the Higgs fields could be 
no other thing as the Top condensate value. [2l|. This 
occurrence could also explain the similarities between the 
properties of the quark mass spectrum and the spectrum 
of superconductivity systems, underlined in the "Demo- 
cratic Symmetry Breaking" analysis [2J. 

In order to evidence the behavior of the potential at 
larger values of X, in Fig^J almost the same plot but 
for an increased range of the variable X is depicted. The 
picture clearly show that, in the framework of the present 
approximation, and for reasonable values of the coupling 
( a = 0.6), there are no terms that control the instabil- 
ity for the generation of the quark condensate parame- 
ter, which under the shown potential will tend to grow 
without limit. Therefore, it becomes clear that the sta- 
bilization of the minimum of the system should come 
from terms higher than the ones considered here. Pre- 
cisely in |l3[ . this behavior was guessed to occur thanks 
to the color interaction between quarks. Therefore, un- 
der the assumption that the technique being used in this 
exploration is well describing the massless QCD, it seems 
that this theory could dynamically develop heavy quark 
masses. This outcome could be another realization of the 
dimensional transmutation effect |22| . A requirement for 
the next corrections to produce helpful results for mod- 
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00 



FIG. 6: The potential plotted as a function of X and g for 
fixed values of M = 333 MeV and /i = 6.8 GeV. It can be 
observed that for each value of X, there is critical coupling 
g below which the potential becomes negative, that is, lower 
than its value at vanishing condensates. This critical coupling 
decreases when X grows. 



elling, is that the stabilizing potential at large X values 
behaves in such a way, that its dependence on M assures 
that the extreme point occurs at low values of M. Then, 
it could be expected that it can be fixed to the observable 
value near 333 MeV by selecting appropriate values for 
the coupling and the scale parameter. Also the value of 
X ~ 175 GeV should be allowed to be fixed. 

The Fig. [3] shows the value of the imaginary part Vi m 
of the potential as a function of M and X. Note that 
the dependence in X is not rapidly growing, a behavior 
that if maintained for large X values and in higher ap- 
proximations will indicate an increasing stability of the 
vacuum being proved, for the interesting region of high 
values of X. The picture is for / = — (|) 2 , as remarked 
before, for the positive value of / the imaginary part of 
the potential vanish. More generally, it can be remarked 
that all the other types of pictures shown in this section, 
after being plotted for the positive value of /, show a very 
similar behavior. Further Fig. [S]show the dependence of 
the potential on the variable X and the gauge coupling 
g. Here the mass parameter M was fixed to 333 MeV 
and again the fi is taken as 6.8 GeV which fixes the min- 
imum in the variable M at X — to be near the value 
of M . It should be recalled that we are considering that 
only one quark is being condensed. Therefore, the light 
quarks which in the present discussion do not develop 
their own condensates, should show the observable value 
of its constituent masses. Since this quantity is fixed by 
the value of M, the graphics selected to be evaluated are 
always chosen to show a minimum near M — 333 MeV 
at X = 0. The picture shows how the potential becomes 
negative (lower than its value at zero condensate state) 
when the coupling increases its strength over an amount 
fixed by the value of the quark condensate X. The greater 
the value of X smaller becomes the critical coupling. The 
Fig0show the dependence of the potential on the gluon 
condensate and the coupling constant for fixed values of 
X = 1 and fi= 6.8 GeV. It can be seen that below cer- 




FIG. 7: The plot of V as a function of M and g for fixed 
values of X = 1 GeV and /x = 6.8 GeV. It shows that there 
is critical value of strong coupling g below which the potential 
at M = grow for non vanishing gluon condensates (M > 0). 

tain critical coupling value near to 2 for all values of M, 
the zero gluon condensate state is stable. Increasing the 
value of X is not destroying this property and the value 
of the critical coupling is simply diminishing for larger X 
values. 



VI. SUMMARY 

The implication of a modified perturbative expansion 
for QCD are further investigated. Firstly, an scheme for 
making well defined the diagrams of the proposed expan- 
sion is introduced. After that, the zero and one loop con- 
tributions to the effective potential are evaluated. Fur- 
ther, the results of the zero and one loop calculation are 
improved by adding the quark condensate dependence of 
a relevant two loop term involving the quark propagators. 
The evaluated potential, in the considered approxima- 
tion, indicates an instability of the masslcss QCD upon 
the generation of quark condensates. At this approxima- 
tion also there is no terms making the potential bounded 
from below. Thus, next corrections should produce such 
terms. Therefore, the results, could be detecting a pos- 
sibility for the identification of a sort of Top Condensate 
model as possible effective action for massless QCD. At 
this point it seems useful to remark, that the source for 
the indicated here effect could not had been yet detected 
through numerical studies, possibly since lattice QCD re- 
sults are still limited in the consideration of the fermion 
determinants. Some questions concerning the gauge in- 
variance implementation and the renormalization proce- 
dure in the scheme are commented as well as some further 
more detailed investigations expected to be considered. 
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APPENDIX A: REGULARIZATION OF THE 
SINGULAR TERMS 

The main technical difficulty for the implementa- 
tion of the modified expansion proposed in the works 
IE El E3 is related with the fact that the addition 
to the standard free propagator term (let us call it below 
the "condensate" propagator) is given by a Dirac's Delta 
function of the momentum. This circumstance, then can 
produce singular diagrams even after the theory is dimen- 
sionally regularized. These singularities are associated to 
the appearance after some loop integrals are performed 
of Dirac Delta functions, or standard Feynman propaga- 
tors evaluated at zero momentum. These factors occurs 
due to the conservation of momentum in each vertex . 
Let us consider a vertex with, let us say n {n = 3, 4 for 
QCD). Then, when n—1 different condensate lines join to 
this vertex, the momentum conservation at it forces the 
value of the momentum at the only resting line to vanish. 
Therefore, if a condensate line is attached to this ending, 
Delta functions evaluated at zero momentum will appear. 
On another hand, when a usual free propagator is con- 
nected to this leg the singularity appearing corresponds 
to 3 , . evaluated at zero momentum. This situation 
should be solved before a full sense could be given to the 
modified expansion. Below, we propose a way of consid- 
ering this problem which clearly should be the subject of 
further investigation for its consistency. Let us consider 
the types separately the two types of singularities 

1. 5(0) singularities 

A direct idea that can come to the mind after consider- 
ing the appearance of these kind of terms is the following. 



As we will employ dimensional regularization, the Delta 
functions in the propagators should be also considered as 
dimensionally regularized forms of the Dirac Delta func- 
tion. However, it has been recognized that similar Delta 
functions evaluated a zero spacial coordinates appearing 
in gravitation theories can be analytically extended to 
continuous dimensions D and moreover, their expression 
after taken in the limit of the real space D — > 4 tends 
to vanish. It can be noticed that this is not a counter- 
intuitive result. This is simply because it is possible to 
impose on the succession of functions defining the Delta 
distribution, the condition of to vanish at the supporting 
point without destroying the possibility that the limit 
of the integral of any continuous function multiplied by 
the elements of the secession, tends to the value of the 
function at the support. Therefore, we follow the same 
procedure here and interpret the Delta functions appear- 
ing as D dimensional ones. Then, it is possible to step 
by step rep roduce the arguments of Capper and Leib- 
brandt [lj| to conclude that these factors should vanish 
after removing the dimensional regularization. Let us do 
it below for the sake of concreteness. 

For the singular $(0) we can write 




This is a singular D dimensional integral in Euclidean 
momentum space, bein g co mpletely similar to a one in 
real space, considered in [l£j ■ Then, it can be also written 
as follows 

r d P D r d P D p 2 
JeWF ~ JeWWp 1 

f°° f dp D p 2 

But, employing the redefinition of the generalized Gaus- 
sian integral for continuous values of the dimension D 
constructed in Ref. [lfj, it is possible to write 



dp 1 



■ exp(— s p 2 ) = 



,\ D exp(~s f&), 
(47t)tt 2 



dp D ~- 

(27T) 



P , 
~D eX P(" 



1 



S p 



(47T) 



D 



/(f) 



exp(-s f(—)), 



and 



14 



where / is the function introduced in 0] for extend- 
ing the generalized Gaussian integral formula for non- 
integral dimension arguments. As it should be, these 
functions vanish for all integral values of D. 

Then, after employing the integral definition of the 
Gamma function 



T(z) 



/ dt i z-1 exp(-f), 
J o 



for the regularized form to the delta function at zero 
momentum follows 



6(0) 



(4tt)* 



2 K 2 1 



■r ( .-f, 



which exactly vanish in the limit D — > 4 + . 

Therefore, we will interpret that the evaluations asso- 
ciated to the modified expansion are done by using the 
above representation for the factors 5(0). Thus, as a con- 
sequence, it will be considered that all the diagrams in 
which only such kind of singularities appear will vanish 
in dimensional regularization. 



b) 



at p = singularities 



For this kind of singular behavior, let us follow the 
physical notion about that the modes at zero momen- 
tum are appropriately described only by the condensate 
propagator. Therefore, it seems reasonable to also reg- 
ularize the dependence of p !i^_ ie at exactly zero value of 
the momentum in a way that vanish at this single point 
in momentum space p = 0. 

Following this idea, let us adopt the following particu- 
lar regularization satisfying this criterion 



1 



ir 



p 2 + ie 



p 2 (p 2 



6 2 ' 







This expression vanish at p — and in the limit S 
leads to the scalar Feynman propagator. 

Let us consider now the general loop expansion after 
assuming the above two prescriptions and the original 
free propagators in the Feynman gauge employed in [lflj 



G a q L(p,m) 



G f q lf2 ( P ,M,S) 



p 2 + ie 



iC 5(jp), 



(Al) 



x ab (p) 



p 2 + ie 

5 ab 

p 2 + ie' 



i S flh C f 5(p), 



and the standard vertices of QCD. Then, whenever the 
6— regularization is employed as described above, it di- 
rectly follows that all the diagrams having a fixed number 



of loops showing both types of singularities will vanish by 
taking the limits in both regulators 



5 
D 



0. 
4. 



Therefore, after taking the limit 6^0, before the 
one D — > 4, the remaining finite diagrams (due to the 
dimensional regularization) can be evaluated using the 
non distorted propagators IjAljl . 

At this point is useful to underline some properties of 
the original diagrammatic expansion based on the above 
propagators, assumed the above explained regularization 
conditions have eliminated all the singular contributions. 
Then it follows: 

1) The appearance of a number m of the condensate 
propagators within a n-loop diagram will eliminate m 
of the n loop integrals associated to this contribution. 
Therefore, the considered diagram will be now an "effec- 
tive" n — m-loop one. 

2) After expressing the condensate parameters C and 
Cf in favor of the ones: 



6g 2 C 

'(27T)*' 



(A2) 



Sf 



9 2 C F 
4tt 4 



Cf, 



which also incorporate a power of order two of the cou- 
pling constant g, it also follows that the n-loop dia- 
grams of the effective expansion, as considered as mul- 
tiple power expansion in the three parameters m 2 , Sf 
and g, also shows the property that, given the number 
of external legs of the diagram, the number of loops is 
fixed by the power of the coupling constant appearing 
in it. This follows directly from the fact that each time 
that a condensate line appears a loop integral is annihi- 
lated and correspondingly the power of g of the diagram 
is reduced by two. Let us consider that p is the power 
of g corresponding to a n loop diagram in which the pa- 
rameters are not redefined. Thus the new power of g of 
this diagram in which m condensate lines and the new 
parameters are introduced will simply be 



P 



p — 2m. 



(A3) 



This property, then seems to allows for a useful reorder- 
ing of the perturbation expansion. To see elements sug- 
gesting it, let us consider a particular n-loop diagram 
in which the change of the parameters l|A2(l have been 
introduced and corresponding line symbols have been in- 
troduced for the standard and the condensate propaga- 
tors separately. Therefore, for any particular standard 
type line in this diagram it seems possible to consider 
the infinite summation of all the zero order in g (tree) 
contributions to the connected propagator (which by con- 
struction have the same number of loops but comes from 
higher loops in the original expansion). This is done by 
considering fixed the other standard lines. Thus, if not 
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blocked by some difficulty associated to the combinato- 
rial and symmetry factors in the diagrams, the performed 
infinite additions seems that can be shown to be possi- 
ble for all the normal lines. These zero order in the new 
expansion propagators are no other things that the ex- 
pressions 0-© written in the first section. This was 
effectively shown in Ref. 01 . Thus, it seems possible to 
demonstrate that the loop expansion can be reordered to 
produce other version of it in which the new tree propa- 
gators will be ©-I©. The investigation of this possibility 
is expected to be considered in further works. 

APPENDIX B: APPENDIX B 

Let us comment here the on the question about the 
implementation of the gauge invariance in the proposed 
scheme. This problem can be divided in two main lines: 
a) The satisfaction of the Slavnov-Ward-Takahashi iden- 
tities given a fixed quantum gauge condition, and b) The 
invariance of the physical quantities under changes in the 
quantum gauge condition, in particular for different val- 
ues of the gauge parameter a [2^ ■ Both issues need for 
additional attention within the modified expansion under 
study, since the new appearing elements, the condensate 
propagators, being neat distributions, make the discus- 
sion more subtle than in the normal situation. This study 
is planned for a next more basic work. However, some 
points of interest can be remarked here below. 

In connection with the satisfaction of the Slavnov- 
Ward-Takahashi identities, already in the work Q, it fol- 
lowed that the one loop correction to the polarization 
operator satisfies the simplest identity, that is, is exactly 
transverse. This is a non trivial result that suggests the 
possibility of its occurrence at higher loops approxima- 
tions. . AlsoinRefs. [Tol fill some general argues and 
particular checks of this property in particular processes 



were given. The study of this problem will be continued. 

On more general grounds, an observation that also in- 
dicates the possibility of implementing the invariance un- 
der the changes of the quantum gauge condition is the 
following one. 

It can be noticed that due to the following identities 
in the sense of the generalized functions: 



the propagators of the modified expansion are also inverse 
kernels of the second variational derivative of the tree 
level action of massless QCD. Therefore, formal steps can 
be done to transform the full generating functional of the 
Green functions of massless QCD as written in the Wick 
expansion representation, to a functional integral repre- 
sentation over the gluon, quark and ghost fields. Hence- 
forth, there is the possibility that some of the changes 
of variables which are employed to show the quantum 
gauge independence of physical quantities could be also 
implemented in the case of our interest, since the func- 
tional integral will only differ in the boundary conditions. 
A particular interesting way of considering the problem 
seems to employ the Yokoyama modification |24| of the 
Nakanishi-Lautrup B field [25| quantization of the inter- 
action free version of massless QCD. In this scheme the 
variation of the gauge parameter a can be implemented 
as a quantum gauge transformation between the field op- 
erators. Thus, the transformation properties of the free 
generating functional associated to the condensate states 
employed in Ref. @, under the gauge parameter mod- 
ifications could be more efficiently investigated. These 
issues requires careful study that we expect to be able of 
perform in next works. 
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